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O ■ Abstract 

I The paper contains the proof that the diffusion ensemble of point wise parti- 

ta ' cles with the intensity depending on the grain of spatial resolution serves as the 

satisfactory approximation of one quantum particle dynamics. 



a 

1 Statement of problem 

The work continues the algorithmic approach to quantum mechanics that has been ini- 
^ ■ tiated in the works [1], [3]. The main aim of algorithmic approach is to build the video 

film reflecting the dynamics of many quantum particles. We believe that the modern 
computers make possible to create such video films that reflect the dynamics of billions 
quantum particles including all known quantum effects. Such video films must show all 
details of process that could be obtained by any sequence of observation fulfilled on the 
considered system. The creation of such picture we call the simulation of quantum dy- 
namics. The simulation can be based on quantum theory only, because it is the single 
theory of micro world. At the same time the simulation of quantum dynamics represents 
the serious difficulty due to the fundamental difficulty of quantum computations. This 
difficulty lies in the formalism of quantum theory: the dimensionality of the space of 
states grows exponentially when the total number of particles grows. 

Though, the computational difficulties arise already for one quantum particle. Let its 
configuration space be divided to N elements. It means that the space of quantum states 
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has the dimensionahty A^. If we apply the matrix algebra in any form for computations 
we force the computer to process all trajectories of the system passing through all basic 
states. In the simplest case it is expressed in the multiplication of unitary matrices of 
the time evolution. The mean value of the module of matrix element is 1/ \/N . If all 
the interference arising in the evolution of the system is constructive we would obtain the 
resulting matrix which coefficients are about 1/{\/N\/N = 1, whereas they must be 
of the order 1/\/N. It means that the bulk of interference is destructive and the huge 
portion of the computational recourse is spent to verify only that there is no particle in the 
considered point. The computational methods of matrix algebra when applied to quantum 
mechanics a priori require the no efficient expenses of the computational resources. Such 
methods factually realize Ryman scheme of integration for Schr?dinger equation that is 
based on the division of the configuration space to finite elements. 

We propose the alternative approach corresponding to Lebeg scheme of integration. 
We take as the basic elements the point wise samples of the considered particle so that 
each sample will represent the whole particle in one point. The total number of samples 
must not be large and its dynamics must be the good approximation of the quantum 
dynamics of the initial particle. Our main aim is to avoid the huge non effective expenses 
of the computational resources featured to matrix algebra. After the developing of this 
technique for one quantum particle we could hope to apply this method for many quantum 
particles as well. 

This method of representation of quantum particle by its samples we call the method 
of collective behavior. It realizes the requirement of strict economy of the computational 
resources, which is the basic principle of algorithmic approach. The requirement of maxi- 
mal economy of the computational resources is not only esthetic. This requirement allows 
to build the models in which decoherence is inbuilt feature of the model but not an axiom 
as in the Copenhagen quantum mechanics. In is shown in [4J how this requirement gives 
us the classical urn scheme for Born distribution of probability for the results of quantum 
measurements. 

Here we give the interpretation of the dynamics of one quantum particle by means 
of collective behavior. The cost that will be paid for the economy of the computational 
resources is the necessity to build the algorithm by means of mechanism of interaction 
between the samples of particle, not by the differential equations. For the method of 
collective behavior it is impossible to build the adequate differential equation. However, 
this situation has the positive sides as well, beyond economical computations. The model 
of quantum dynamics becomes nearer to classical than in the standard approach that 
makes possible its visual representation. 

The proposed approach is the direct generalization of diffusion Monte Carlo method 
to the case of the time dependent solution of Schr?dinger equation. The known fact 
that DMC gives the most exact approximation of stationary wave functions among all 
computational methods inspires optimism in the practical application of the method of 
collective behavior for more complex problems. 
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2 Dynamical diffusion swarm 



Here we define the main instrument of quantum simulation: dynamical diffusion swarm. 
This object generalizes two well known notions: the ensemble of point wise particles from 
DMC method, and the ensemble of particles with the interaction induced by some classical 
Hamiltonian H{r,p). Particles from DMC have no speeds and they are designed for the 
computation of stationary states for which they give the best approximation. The density 
of particles p for the ensembles with classical Hamiltonian depends on the coordinate r 
and on the impulse p; it obeys Liouvill equation 

The behavior determined by this equation cannot simulate quantum evolution with the 
admissible accuracy because it does not give principle quantum phenomena like Rabi 
oscillation or quantum spectra. Hence, for the simulation of quantum mechanics we must 
admit some elements of behavior of the samples which do not follow from the classical 
physics. 

The next analog of the dynamical diffusion swarm is Calder Legett model for the 
partial decoherence of quantum particle in which the particle is considered as interacting 
with the bath of harmonic oscillators. This interaction gives some random speed to the 
particle. However, such a model is based on standard formalism whereas the dynamical 
diffusion swarm is designed to replace this formalism. 

Why the dynamical diffusion swarm is better than the immediate solution of Shroedinger 
equation ? In this solution we factually use Riemann schem of integration. We niTist per- 
form computations of the wave function on the whole configuration space independently 
on the degree of constructiveness of the interference. Here on the main part of the space 
where the interference is destructive and the wave function factually equals zero we are 
forced to spend the computational recourse only to verify this. The dynamical diffusion 
swarm, in contrast, realizes the more general Lebeg scheme og integration. In this case the 
diffusion dynamics results in that the samples will concentrate in the areas of construc- 
tive interference themselves and we avoid the non effective expenses of the computational 
resources. This is the fundamental advantage of the diffusion dynamics. We will see that 
the cost for this is the non uniform dependence of diffusion rate on the grain of the length 
5x in contrast to the standard diffusion where the rate is uniform. 

We proceed with the definitions. We call the swarm the finite set S consisting n 
identical point wise particles each of which s E S has its own coordinates and impulse 
x{s), p{s) e R^. In the method of collective behavior one quantum particle of the mass 
M and charge Q is represented by the swarm S each member s E S oi which has the mass 
m — M/n and the charge q = Q/n. The elements of this swarm are called the samples 
of this quantum particle. We suppose that the total number of samples n is so large that 
the swarm can be used as the approximation of the continuous media. E.g., if we need 
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use the smaller and smaller spatial grain some samples will always occur in each spatial 
cell. But the dispersion of speeds will grow when the grain decreases, and we will have 
the separate swarm for each spatial grain 5x. 

The choice of spatial grain is closely connected with by the definition what object 
must be considered as quantum particle. Thiis definition in turn depends on the concrete 
problem and quantum particles are not necessary elementary in the sense of theoretical 
physics. The definition what must be treated as a particle presumes the choice of the 
typical length AX and the time AT, so that the size of particle is much lesser than AX, 
e.g. it can be treated as point wise, and the time interval AT is not less than the typical 
time of the processes we are interested in. Let us agree that the typical mean speeds of the 
considered shifts are much lesser than some limit speed of all movements c. For example, 
an atom can be treated as a point wise particle in he processes with /S.X > lO^^m and 
AT > 10~^°s. If we decrease the value of typical lengths and times then to obtain the 
right picture we must consider the different set of elementary particles, for example, the 
separate nucleus and electrons inside of atom. If we fix AX and AT, then to obtain the 
dynamical picture we must define the smaller segments 5x, 6t, which will represent the 
elementary steps of the video film and which, however, must be much greater than the 
typical lengths and times AX, AT of the more fundamental processes then the considered 
one (the gap between the different fundamental processes can be about 10"^'', that always 
allows to make this separation). Also in the same process with the fixed energy the times 
and lengths depends on masses. The separation of particles by their masses makes possible 
to consider for the bulk of processes in electrodynamics only electrons because the typical 
distances of flight of protons are to 1800 times smaller. We can then treat the chosen 
values AX, AT as the size of imaginary screen and the length of video film, and 6x, 
6t as the grain of spatial resolution of screen and the time of showing of one card in 
the film. We choose 6x and 6t maximal so that our film will be informative. After this 
choice the conclusion can be done about what particle should be treated as quantum and 
what as classical. For this the typical values of their action a = M(AX)^/AT should be 
compared with Plank constant h. If a < h then the particle should be treated as quantum, 
otherwise as classical. In the method of collective behavior the passage from classical to 
quantum type of consideration means simply the change of swarm size, e.g. does not 
mean the different type of dynamics. Due to the above mentioned reserve in the choice 
of resolution in the process of film preparation we then can further decrease the values 5x 
and St for the forming of right picture , for example, dividing these intervals to smaller 
parts and obtain the better approximation to the solution of Shroedinger equation. We 
assume that the space is divided to the equal cubes with the side 6x, and the time is 
divided to the equal intervals of the longitude St. 

We introduce the value c, which is the single nonzero speed of movement of the samples 
in the swarm. The intervals of time and distances will be always chosen so that Sx S> cSt. 
It guarantees that in each step of the evolution the values of magnitudes obtained as 
the mean values on cubes with the side Sx will change small that is necessary for the 
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asymptotical approximation. 

The density of swarm in the point x is determined by the expression 



p(r,t) 



N{r,t) 



(1) 



where iV(r, t) denotes the total number of samples occurring in the moment t in the same 
cube with the point r. For the comparison with Shroedinger equation in this definition 



with the densities pn with the growing n instead of one swarm. Further we assume that 
the value 5x is fixed. We write p{x) = |\E'(a;)p instead of 



where the convergence is uniform without special mentioning. This sequence of the swarms 
realizing the approximation of the exact solution of Shroedinger equation is called the 
admissible approximation of quantum evolution. 

Our aim is to define the behavior of the samples in the swarm which gives the admis- 
sible approximation of quantum evolution. 

The main requirements to the simulation of quantum dynamics through the collective 
behavior are the following. 

• Quantum dynamics is simulated by the dynamics of the swarm of samples so that 
in each time instant t the quantum probability equals the density of the swarm. 



in each point x of the configuration space. 

• Each sample of the swarm has its own history, e.g. it preserves its individual num- 
ber in course of the whole simulating process. The types of the samples exactly 
correspond to the types of real physical particles. 

• The behavior of each sample is completely determined by its own state and the state 
of all samples in its close vicinity. 

The swarm satisfying these conditions is called the quantum swarm for one quantum 
particle. 

We define the behavior of the samples such that these conditions are satisfied. For this 
it is sufficient to show that for the solution t) of Shroedinger equation it is possible to 
move the samples only locally, e.g. to the small distance for the insurance of the equation 
([3]) in each time instant. Such a movement, of course, will be a priori non natural in the 
dynamical sense, but we will show how it can be done by means of the dynamical diffusion 
mechanism. 



we should converge 6x 



0, that means the consideration of the sequence of the swarms 




(2) 



^{x,t)\^ = pix,t) 



(3) 
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We note that the second rule means that we refuse from the using of complex numbers 
in the description of quantum mechanics. Also the locality of all interactions allows 
including QED to our model. The behavior of samples is the rule determining the change 
of its internal states (the type, impulse, momentum of impulse) and the spatial position 
(spatial shift). In view of the above mentioned the behavior cannot be determined by the 
classical physics. 

We define the quasi classical behavior of the samples called the dynamical diffusion 
mechanism. The swarm of samples with such behavior satisfies these conditions. 

Let us agree that each sample in each time instant can either to stay in place or to 
move along one of the coordinate axes OX, OY, OZ with the speed c. 

We call the reaction of change the sequence of the following operations on the swarm: 
the choice of pair a, (3 of the samples located not farer than 5x from each other, which 
speeds are mutually opposite: v{a) = — f (/3) and either simultaneous replacement of their 
speeds to zero (if they are nonzero) or acquiring them mutually opposite speeds of the 
module c oriented along one of the coordinate axe. The axe is always chosen randomly 
from the uniform probability distribution. 

The reaction of change does not change neither summed impulse of the swarm, nor 
summed momentum of impulse if Sx is sufficiently small. By N{r) and Ns{r) we denote 
the sets of all samples in the cube with the point r and the set of samples from this cube 
possessing zero speeds correspondingly; by (r) , (r) , (r) we denote the sets of 
samples from the cube r, moving along the corresponding axe in the positive direction, 
and by analogous symbols with the sign — the corresponding sets but for the negative 
direction of movement. By \g\ we denote the total number of samples in the set g. Let 
us agree to denote the total numbers of samples in a set by the same symbol as this 
set but with the replacement to n. We call r— stationary each subset S C n(r) 
consisting of the samples with nonzero speeds for which Ylaes'^i'^) ~ ^ ^"^^ ^ 
maximal on including with this property. The total number 15*1 of the samples of r— 
stationary set (which does not depend on its choice) is denoted by s(r). Let d > he a. 
chosen constant such that the coefficient of the diffusion is proportional to d, V{r) is a 
scalar field proportional to the external potential energy with the constant coefficient of 
proportionality, grad V{r) = (yx{r),Vy{r),Vz{r)). 

We also agree to consider only non relativistic swarms, e.g. such that ns{r)/n(r) is 
close to 1 for all r. It means that the bulk of samples in each cube have zero speed. This 
requirement is incompatible with the point wise approximation by the swarms [2] of the 
exact wave functions for the external potential of Coulomb form 1/r because the mean 
speed of samples in the vicinity of zero point for such potentials must converge to infinity. 
For the asymptotic convergence [2] we would have to assume that c can be chosen as large 
as needed for every next swarm number n. In the reality c cannot exceed the speed of 
light that establishes the natural limitation to the accuracy of the swarm approximation 
of the solutions of Shroedinger equation. 

The dynamical diffusion mechanism of evolution is the following sequence of the op- 
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eration on the swarm: 

• 1) The sequence of random reactions of change with the uniform distribution of 
probabihty, which gives the distribution of the speeds with the property s(r) / ns{r) = 
d for each point r. If n(r) is small, this equality must be satisfied with the maximal 
accuracy (see agreement about the accuracy from above). 

• 2) The acquiring of speeds to some samples from Ns{r), randomly chosen from 
the uniform distribution so that for each axe the signs of newly acquired speeds 
along this axe are the same, and if f„(r) is the summed speed acquired to the 
samples from r-th cube along the axe u, u = x, y, z, then for all such u the equation 
Vu{r)'m = —Vu{r) is fulfilled with the maximal possible accuracy. 

• 3) The change of coordinates r{a) of each sample corresponding to the Galileo law: 

rnew{a) = r(a) + v{a)At. 

• 4) Recalculation of V{r) accordingly to the new positions of samples. 

We do not concretize the method of recalculation of the potential energy. It may be 
done by Coulomb formula or by the diffusion mechanism as was proposed in the work . 

The dynamical diffusion swarm cannot be represented as an ensemble of point wise 
particles with the classical interaction. The item 1) says about two things: 

• there is the force with the random direction which acts to the samples (compare 
with [1]), and 

• the samples acquire the mean value of speed within the accuracy determined by 6x 
(the lesser 6x is the more accurate mean value is taken). 

For each time instant t if Ax is sufficiently small then the density of swarm p{r, t) for 
any point r will not depend on the orientation of the coordinate axes. Indeed, let 6i be 
such that c5t <^ 5ix -C and let v{r) denote the average speed of the samples in the 
point r, found by the averaging on the samples with the coordinates ri : ||r — ri|| < 5iX. 
The total number of the samples occurred in the unit of time from the vicinity of the 
point Ti to the vicinity of the close point r2 will be then proportional to the dot product 
^^('"i)('"2 — '"i)/||'"2 — which does not depend on the orientation of the coordinate axes. 

A state of the dynamical diffusion swarm is determined by the coordinates and speeds 
of all its samples. This is its principal difference from the ensemble of DMC method where 
there are no speeds of the samples. 
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3 Differential equations for the dynamical diffusion 
swarm 

There is no system of differential equations on the density of the diffusion swarm and its 
average speed which is equivalent to Shroedinger equation. 

Nevertheless, there is the sequence of such systems which realizes the admissible 
asymptotic approximation of the solution of Shroedinger equation. Each of the system 
from this sequence depends on the fixed elementary length 5x. For example, the inten- 
siveness of the diffusion process will beb proportional to {5x)~^ . This does not allow to 
launch 5x to zero as it is always done in the mathematical analysis when it is applied to 
the processes of classical physics. The value of grain 5x must be chosen such that the 
approximation of the density p = |\E'p of the wave function by the density of the diffusion 
swarm within {5xY is satisfactory for the considered process. Only after the fixation of 
5x it is possible to build the diffusion swarm of the corresponding intensiveness and the 
differential equations approximating its dynamics which will be equivalent to Shroedinger 
equation. 

A state of the dynamical diffusion swarm is determined by the pair of functions 

p(t,f), p(t,f), (4) 

where p is the scalar function of density of the samples, p{f) is the vector function resulting 
impulse of the samples in the point f defined as lim P{r,dx)/{dxY, where P is the 

dx >0 

summed impulse of the samples occurred in the cube with r with the side dx. Here 
we assumed that dx can be done sufficiently less than the grain 6x, determining the 
coefficients of the equation on p and p. 

The dependence of equations on the grain Sx will be revealed as follows. The summed 
impulse p{t, f) varies slowly when f changes on values larger than Sx. But its derivative 
^ will be very large: of the order l/{6xY, and will vary rapidly as well. E.g., the graph 
of the function p{t, f) is sufficiently smooth if we look at it with the large grain Sx, but if 
we raise the resolution by decreasing the grain 6x, we see that the graph looks like a saw 
with acute teeth. The more is resolution l/6x the more is sharpening of the teeth, it is 
limited by the limit of speed c (compare with [Ij). This important for further condition 
we call the non relativistic approximation and will write it as f <^ c. 

In view of isotropy of the diffusion process the change of density p(r, t) in the time and 
its second derivative can be found by the integration on the surface of the sphere S{r) of 
radius Sx by formulas 

= / p{r,t)n{fi)ds{ri), 

%i) J'r-^n{n)ds{n). 

S{r) 

These formulas are right for any mechanism of changing of the speed of samples. 
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Now we derive the law of changing of the resulting impulse of the swarm in the 
small sphere with center in the point f, which results from the moveme of the samples 
along the vector a of normal to the surface of sphere of the unit length. Three magnitudes 
make the deposits to the change of the resulting impulse: 

• Penetration of the samples which have acquired the speed in the reaction of change 
through the small element of the surface (diffusion). 

• Penetration of the samples which have acquired the speed from the action of external 
potential. 

• Penetration of the samples which have preserved their speed (by inertia). 

It follows from the definition of the diffusion process that these deposits equal corre- 
spondingly — / grad pa, —Kp grad Va and gppd, where I,K,,g is the intensities of the 
corresponding processes. The choice of units system allows us to make g = 1. The depen- 
dence of the grain of spatial resolution needed for approximation of Shroedinger equation 
has the form: 

J - , - ^ . (6) 

2m'^(5xy^ mSx 

In view of the non relativistic supposition we can omit the last summand which is 
sufficiently smaller than the first two for the small 6x. We then obtain the following 
approximate formula. 

dp 

— ~ — / grad p — Kp grad V. (7) 
The resulting equation on the density of the diffusion swarm has thus the form: 
d'^ p(r) f 

= — I grad p — K,p grad V)n{r') dS{r'), (8) 

5(r) 

where the coefficients /, k can be found by [61 

We prove that the quantum swarm satisfies [HI which means the admissible approxi- 
mation of the quantum dynamics by the evolution of the swarm. 

The method of collective behavior permits to give the simple algorithm for the compu- 
tation of the energy, impulse and momentum of impulse of quantum particle represented 
as the swarm of samples. Here accordingly to the quantum rules for the finding of any 
magnitude we have to averaging on the values of the reciprocal magnitude. To find the 
impulse we must fix some value of the time interval At and fulfill averaging of impulses 
on all samples on any values of the distances they overcome. This receipt in our notations 
gives the vector 



n^{r)+ - na;{r) nyir^ - ny{r) n^{ry - n^{r) 
cm — — — -- — ^^,cm-^^^ — , , , cm — — — , , ) (9) 

n[r) n[r) n[r) 
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which equals to the average impulse of all samples in the swarm found by averaging on the 
passes in the fixed time interval. The analogous calculation of the momentum of impulse 
or the potential energy gives the average momentum or the average potential energy. 

When calculating the average kinetic energy we have to fix the pass 5x and to take the 
average energy on all instant of the time t, because the time is the reciprocal magnitude 
for the energy. It means that we must sum only the energies of the samples which are 
moving in the considered time instant, e.g., for any cube for the speed v — cux/n along 
the axe x the expected total number of such samples is nv/c, and their fraction in the 
total number of samples is rix/n. Here we use the non relativistic assumption that the 
total fraction of the moving samples is small. The total kinetic energy found by this rule 
is 

m(?nxnx/n + m(?nyny/n-\- mc^nzUz/n 

which coincides to the kinetic energy found by the conventional formula Af u^„g„„/2 = 
{{crix/nY + {cny/nY + {cnz/nY)nm/2. The laws of conservations for the impulse, mo- 
mentum of impulse, energy of the swarm then follow from the classical laws of conservation 
and the non relativistic assumption. In the next paragraph we prove that the diffusion 
swarm dynamics can give the admissible approximation of the quantum one. Using Eren- 
fest theorems and the laws of conservation for quantum dynamics we conclude that our 
method of calculation of these magnitudes A gives their quantum mean values found by 
the formula {A) = / '^*{r)A'^{r)dr. 

4 About the diffusion swarm with non uniform in- 
tensity 

The intensiveness of the diffusion equals the coefficient of the Laplace operator in the 
diffusion equation. The intensiveness determines the total number of samples passing 
through the unit of square in the unit of time. To simulate quantum dynamics we need 
the diffusion swarm with the non uniform intensiveness. It means that the intensiveness 
of diffusion depends on the chosen grain of the spatial resolution 5x. In this section we 
discuss how non imiform intensiveness of diffusion can be obtained in the same swarm. 
The concrete mechanisms discussed here will not be used further; they can be interesting 
only for the programming realization of the dynamical diffusion method. 

We at first consider the case when the external potential is constant grad V — 
The dynamical diffusion process with the non uniform intensiveness can be insured by 
the special mechanism which we call threads. We illustrate the method on the following 
example. We suppose that all samples move not in all the space but along one closed 
trajectory (thread) which is determined by the smooth embedment of the circle to the 
space: 7 : — > R^. We suppose that the change of speeds happens only along 
this trajectory so that the samples remain in this thread in each time instant. This is 
equivalent to the imposing of the holonomic tie to the samples. We then suppose that the 
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linear density as well as the module of speed of the samples is almost equal in all point of 
the trajectory. We consider the cube containing one point of this trajectory. The flow of 
samples through its border will not depend on its size 5x, because the thread is only one. 
The intensiveness of this process is thus proportional because the quantity of 

samples penetrating in unit of time into the cube with the side 5x does not depend on 5x, 
and the density is obtained by the division of the quantity of samples to the volume. This 
example is not very good because many areas of the space remain without the samples at 
all. 

We consider the next example. Let the space be divided to the cubes which are 
grouped by the layers 1,2,.. .. For each j = 1, 2, . . . the cubes of the layer j + 1 consist 
of 8 cubes of the layer j, and their side, correspondingly, is twice large. For each j the 
change of samples between the neighbor cubes of the layer j, occurring to the same cube 
of the layer j + 1, goes only through the narrow channel with small capacity independent 
of j. The quantity of samples moving between the cubes of any fixed layer will not then 
depend on the number of this layer. It can be guaranteed by the appropriate choice of 
the pair for the change of impulses. Such a mechanism gives us the required intensiveness 
of the diffusion proportional to l/(5a;)^, in view of the definition of density 

Now we consider the case of varying external potential. At each step of the evolution 
for the samples acquiring or loosing their speeds in the change we will use the rule from 
above, which insures the intensiveness proportional to l/((5x)^. The samples acquiring 
their speed from the action of the external potential will move as usual, independently 
of layers. We then obtain the formula This space design show how in principle the 
necessary non uniform intensiveness can be obtained in on swarm. 

Practical realization of these methods presumes that we trace for the spatial location 
of the separated parts of the swarm, which means the refusal from the uniformity of the 
space and passage to the fractal space. The space with fractal dimensionality arises if 
we use the non uniform grid for the method of finite elements. This is why the proposed 
computational receipt says that we must fix the grain 5x of the linear resolution such 
that the corresponding approximation of the wave function is satisfactory for our aims, 
and then consider the diffusion swarm with the intensiveness found by [61 If we are not 
satisfied with the obtained dynamical picture, we must choose the new value of 6x and 
repeat all the work. 

5 Equivalence of quantum and dynamical diffusion 
swarms 

Here we show that the sequence of diffusion swarms represents the admissible approxima- 
tion of quantum evolution. We have defined the quantum swarm as the swarm satisfying 
[3] and which evolution can be represented as the local movements of the samples. 

At first we determine that the quantum swarm exists, e.g., that the equation [3] can be 
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really reached by only local shifts of the samples. Then we prove that the mechanism of 
the movements of samples coincides with the diffusion that gives the main result. 
Let us consider the quantum swarm. We start with Shroedinger equation 



ih 



dt 



2M 



A*(r,i) + V;<,t(r,i)*(r,i) 



(10) 



which can be rewritten as 




_|_AvI/j(r) + ^M/^(r), 



(11) 



for the real and imaginary parts ^'^ , of the wave function ^. We focus on the evolution 
of the density of quantum swarm, which is the function 



Fixing the value of 6x we apply for the approximation of the second derivative the 
difference scheme of the form 



for each time instant, where the wave function is supposed to satisfy all the sufficient 
conditions for such approximation. Since the addition of any constant to the potential 
energy Vpot docs not influence to the quantum evolution of the density, we can consider 
instead of Vpot the other potential V = Vpot + o:, where a = — ,„^'^^'^^2 , that leads to the 
disappearing of the summand 2\l'(a;) in the difference schemes for the second derivative on 
X, y, z (from which the coefficient 3 arises) after its substitution to Shroedinger equation. 
We introduce the simplifying coefficient 



Since we yet do not know the mechanism of moving of the samples in quantum swarm, 
we suppose that we simply either take off some quantity of the samples from any cube 
or put it there from some storage. We split the evolution of the quantum swarm in the 
time to so small segments of the longitude St, that in each of which the samples move 
only in the framework of two neighbor cubes. If we prove that the evolution of quantum 
swarm on such segment is insured by the diffusion mechanism, it will be right for the 
whole evolution as well. That is our supposition does not limit the generahty. We also 
agree that these cubes differ from one another by the shift to 6x along the axe x, that 
does not limit the generality as well. Let the centers of these cubes be x and Xi — x + 5x. 
In these suppositions the summand ^(x — Sx) in the difference scheme disappears as well 



d^'^ix) _ *(x + Sx) + ^(,T - Sx) - 2-^{x) 
dx'^ (Sx)"^ 



h 1 



7 = 



2M (5x)2' 
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and on the considered small time segment the evolution of quantum swarm is determined 
by the system of equations: 



and by the analogous system obtained by the substitution of x in place of xi and vise 
versa. 

For such a segment we thus have 

^ = 2^*(x)(7^'-(a;i) - +2^'-{x){--f^'{xi) - V{x)^'{x)) = 

It results in that the outcome of the samples in one cube equals their income to the other. 
The evolution of quantum swarm then satisfies the condition of locality. In the order 
to compare these evolution with the diffusion we now find the second derivative of the 
quantum density to time: 

= 27[(7*''(a;i) - V{x)^''{x))^'-{xi) + ^%x){--f^'{x) + V{xi)^'{xi))- 
(-7^^(xi) + V{x)<i/'{x))^'lxi) - ^'■(x)(7^''(x) - V{xi)'^^{xi))] = 
272(^'-(xi))2 - 27l/(x)*"(x)^' (xi) - 2-f\^'{x)y+ 

27r(xi)^*(x)^^(xi) + 272(^*(xi))2 - 27l/(x)^^(x)^*(xi)- (14) 
2'-f'^{W{x)f + 2-fV{xi)W'{x)W{xi) = 

272((^'-(Xi))2 + - + + 

27[(V^(xi) - V{x)){{^^{x)f + {"^'{xyf) + o{6x)], 

where o{Sx) = (^'■(x)^^(xi) + - ((^'^(x))^ + {<iJ'{x))'^)){V{xi) - V{x)). Now 

we compare it with the expression for the second derivative of the density of diffusion 
swarm found in the previous section, taking into account that in our case the change 
of samples goes between two neighbor cubes along the axe x only. Comparing with [8] 
in view of [6l we conclude that the second derivative of the density of quantum swarm 
asymptotically converge to the second derivative of diffusion swarm. 

If we choose instead of the initial state the state where its density has Gauss form 
which is the ground state of harmonic oscillator, then for the corresponding value of 
energy V = a(x^ + + -z^) we have dp/ dt = in the initial instant for any point of space. 
It is proved that the second derivative of the quantum swarm density and of the diffusion 
swarm density are the same, hence the diffusion swarm will be good approximation for the 
quantum density on some interval AT. Switching on slowly some potential we will have 
also the approximation of any quantum evolution in the limit of swarms for the unlimited 
increasing of n. 

The swarm approximating quantum dynamics of one particle depends on the choice of 
6x. After the fixation of 6x we obtain for unlimitedly decreasing 6t the approximation of 
the wave function within 6x. The intensiveness of the diffusion will be determined by 6x, 
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it will be yy^3^(^3.)3 • If "we want to decrease the grain 5x we must allow the more quantity 
of moving samples in the unit volume. It is necessary due to the uncertainty principle 
for length and impulse: the dispersion of speeds of the samples will grow if 5x decreases. 
In any case for the obtaining of the dynamical picture one has to fix the grain of spatial 
resolution 5x. 

If the total number of samples n is limited we obtain the model of quantum dynamics 
with decoherence of the inbuilt type. This model can be extended to the multi particle 
case where it can serve as the approximation of quantum dynamics in the standard Hilbert 
formalism (see ([!]). The robustness of this scheme for the numerical computations follows 
from that it gives Born rule for the quantum probability which thus turns to be inbuilt 
into the algorithmic formalism, in contrast to Copenhagen formalism where this rule is 
postulated. 

6 Restoration of wave function from dynamical dif- 
fusion swarm 

We have solved the problem of approximation of the dynamics of density for one quantum 
particle by the special diffusion process with the non uniform intensity. A state of the 
dynamical diffusion swarm is determined by a pair HI Such a pair does not contain 
the notion of complex numbers which induced the famous quantum interference in the 
standard formalism. Furthermore, the diffusion swarm gives no beautiful differential 
equations of Shroedinger type for p and p. It radically differs from the classical processes 
(for example, heat transform or oscillations) because its intensity depends on the chosen 
grain of spatial resolution. We agree to these for the sake of the main: the economy of 
the computational resources required for the description of quantum dynamics. 

Now, to finalize the picture we have to solve the inverse problem: to show how to 
restore the wave function \1/ from the given state of the dynamical diffusion swarm |H to 
do this we turn to the equality [131 ^-iid substitute to it the expression of the wave function 
through the density: \l/(r) = a/ p{r) exp{i(f){r)). We must find the phase 0(r) of the wave 
function. Since only relative phase has the physical sense we can fix some point r and 
consider the phase of the other point ri relatively to r. If ri is close to r, the equation [13] 
gives us 



(ri) = arcsin k{6x) 



2 p{r-ri] 



v/p(r)p(ri) 

which results in the following formula for the finding of the relative phase: 

0(ri) = j k{5xfv (15) 

7 

where the contour 7 goes from r to ri. This definition depends explicitly on the choice of 
contour 7, hence we have to prove its correctness, e.g., its independence of the contour 7. 
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Since the phase is determined within an integer multipher of 27r, the different choices of 
the contour can resuh at most in the change of the phase on such a number that takes 
place in case of excited states of electron in hydrogen atom with the nonzero momentum 
(for example, 3d). We show that the integration of speed v of the swarm on the closed 
contour preserves its value in the time with the more exactness the less 5x is. It involves 
that if in the initial time instant the definition f[T^ was correct, it then remains correct 
for the next time instants. 

We now consider the derivative of the integral of the speed for the closed contour 7c. 
Applying the formula ([8]) and accounting dp/dt = p dv/dt, we obtain 



J {6xyv dl = - J I{5xY^^—^ + k{5xY grad V. (16) 



7c 



The first summand gives zero after the integration on the closed contour, because it is 
grad In p, the second summand gives zero by analogous reason. 

Now it is sufficient to verify the correctness of the definition (fT5|) in the initial time 
instant that can be done immediately for any task. If the wave function of initial state can 
be obtained from the ground state in the Coulomb field where w = 0, then the correctness 
follows from the proved because there is no any phase shift here to 27r/c. If for the obtaining 
of the initial state in the considered problem we must start from some excited state with 
the shift of the phase, the correctness should be checked for this state at first. 



7 Conclusion 

We have given the economical method of the simulation of one quantum particle dynamics. 
This method allows the creation of the video films showing the evolution of its state in 
the time. Our method is based on the mechanism of diffusion of the samples of swarm 
consisting of the point wise particles which obey the laws of classical mechanics with only 
one additional property: they change by impulses. The change of impulses creates the 
mechanism of diffusion for this swarm. The intensity of diffusion depends on the chosen 
grain of spatial resolution: the less the grain is the more is the intensity. This is why such 
a swarm cannot be described by a beautiful system of differential equations. Its main 
advantage over the methods of matrix algebra is the strict economy of the computational 
resources required for the simulation and the similarity to the classical description of the 
dynamics that permits to build the visual pictures of the evolution. 
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